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translates  the  indepth  2-dimensional  temperature  field  and  PGPLOT  which 
plots  the  resultant  finite  element  mesh  and  temperature  contours.  The 
POESSY  code  is  designed  to  accept  the  output  data  from  a thermochemical 
ablation,  shape  change,  in-depth  temperature  response  computer  code  and 
prepare  it  for  input  to  a thermostructural  analysis  computer  code  with 
a minimum  amount  of  effort  by  the  user. 

The  use  of  the  POESSY  computer  program  is  demonstrated  by  sample 
problems  solved  for  a typical  plug  and  shell  nose  tip. 
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The  POESSY  computer  program  was  developed  to  help  fill  a gap  in 
existing  reentry  vehicle  nose  tip  analysis  capability.  Programs  exist 

I 

to  analyze  the  ablating  shape  and  indepth  temperature  response,  such  as 
PAGAN  (Reference  4)  and  to  analyze  the  thermostructural  response  such 
as  DOASIS  (References  1,  2,  and  3).  However,  no  preprocessor  programs 
exist  to  automatically  translate  the  output  from  a program  such  as  PAGAN 
into  an  appropriate  form  for  input  to  a program  such  as  DOASIS.  This 
[ translation  task  is  not  trivial  since  different  forms  of  problem  definition 

are  used  in  these  two  analysis  tools.  Consequently,  POESSY  was  developed 
to  perform  this  translation  process  automatically,  with  a minimum  of 
required  user  input  information. 

The  POESSY  computer  program  consists  of  four  basic  subprograms  which 
perform  different  distinct  steps  in  the  translation  process.  They  are 
PGMESH  which  automatically  generates  the  structural  finite  element  mesh, 
PGPRES  which  translates  the  surface  pressure  and  temperature  information, 
PGTEMP  which  translates  the  indepth  temperature  distribution  and  PGPLOT 
which  plots  the  resultant  finite  element  mesh  and  isotherm  contours. 

This  report  describes  the  theory  behind  these  subprograms,  the  form  and 
type  of  input  data  necessary  to  run  POESSY  (User's  Manual)  and  several 
sample  problems  illustrating  the  use  of  the  POESSY  computer  program. 
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INTRODUCTION 

The  POESSY  computer  program  consists  of  four  subprograms  which 
perform  the  mesh  generation  (PGMESH),  the  surface  pressure  and  temperature 
translation  (PGPRES) , the  indepth  temperature  translation  (PGTEMP)  and 
the  plotting  of  the  mesh  and  isotherm  contours  (PGPLOT) . All  of  these 
programs  except  PGMESH  were  developed  by  simply  modifying  the  existing 
DOASIS  pre-  and  post-processor  computer  programs  PRSINT  (for  PGPRES), 

TEMINT  (for  PGTEMP)  and  CONTUR  (for  PGPLOT) . The  complete  DOASIS 
family  of  computer  programs  is  described  in  References  1,  2,  and  3, 
with  PRSINT,  TEMINT  and  CONTUR  being  described  in  Reference  2.  Since  the 
POESSY  subprograms  are  so  closely  tied  to  these  pre-  and  post-processor 
programs,  the  reader  should  consult  these  references  for  relevant  and 
additional  theory  and  programming  techniques,  since  this  report  uses 
the  information  contained  in  these  references  to  help  describe  the 
theory  behind  the  POESSY  subprograms. 

The  automatic  nose  tip  mesh  generator  PGMESH  is  described  in  Section 
II.  It  was  developed  for  the  specific  purpose  of  generating  nose  tip 
finite  element  meshes  and  therefore  uses  techniques  not  contained  in  the 
DOASIS  general  purpose  mesh  generator  MESHGN  described  in  Reference  2. 
However,  the  overall  philosophy  of  using  an  associated  indical  I-J 
grid  with  the  coordinate  R-Z  grid  is  used  in  PGMESH  and  therefore  Sections 
1 and  2 of  Reference  2 should  be  consulted  for  a better  understanding  of 
this  concept.  Section  II  of  this  report  basically  describes  the  technique 
used  to  automatically  generate  the  mesh  from  the  given  PAGAN  code  input  data. 

The  automatic  nose  tip  pressure  translator  PGPRES,  temperature 
translator  PGTEMP  and  plotter  PGPLOT  are  described  in  Sections  HI#  IV,  and 
V respectively.  These  subprograms  are  direct  modifications  of  the  programs 
PRSINT,  TEMINT  and  CONTUR  described  in  Reference  2,  respectively,  and 
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therefore  will  referenced  directly  when  describing  PGPRES,  PGTEMP 
and  PGPLOT. 


A complete  description  of  the  POESSY  computer  program  is  presented 
in  Section  VI,  followed  by  a description  on  how  to  modify  the  program 
size  in  Section  VII.  The  POESSY  User's  Manual  is  presented  in  Section 
VIII  followed  by  two  sample  problems  given  in  (Appendix  A)  which  illustrate 
the  use  of  the  POESSY  code.  Appendix  B contains  a description  of  plotting 
software  used  by  the  POESSY  computer  program. 


SECTION  II 


PGMESH  AUTOMATIC  NOSE  TIP  MESH  GENERATOR 
The  set  of  subroutines  contained  in  the  segment  PGMESH  are  designed 
to  automatically  generate  a structural  finite  element  mesh  when  given 
only  the  outside  boundary  contour  of  the  solid  body.  This  outside 
boundary  definition  is  all  that  is  needed,  provided  some  criterion 
is  choosen  as  to  how  to  subdivide  the  interior  into  finite  elements. 

The  criterion  choosen  for  POESSY  is  to  use  a "regular  rectangular 
pattern"  of  finite  elements,  with  transition  elements,  that  is, 
arbitrary  triangles  and  quadrilaterals,  existing  along  the  boundary 
contour.  This  concept  of  a regular  rectangular  pattern  of  elements 
is  graphically  illustrated  for  a plug  and  shell  nose  tip  in  Figures  1 
and  2 respectively. 

When  viewing  these  figures,  one  will  notice  that  the  origin 
of  the  rectangular  pattern  is  chosen  to  coincide  with  the  axis  of 
revolution  of  the  nose  tip  for  the  radial  (R)  coordinate  and  with  the 
backside  surface  for  the  axial  (Z)  coordinate.  The  main  reason 
basically  for  this  choice  of  pattern  and  origin  is  that  for  receding 
nose  tips  the  outside  front  surface  is  the  only  surface  which  changes 
as  recession  takes  place.  Consequently,  the  location  of  elements  when 
referenced  from  the  backside  will  remain  fixed  provided  the  size  and 
shape  of  the  basic  rectangular  element  remains  fixed  for  meshes  of 
the  nose  tip  at  different  points  in  time  except  for  elements  on  the 
surface  which  are  effected  by  the  recession  of  the  surface.  Therefore, 
one  can  trace  the  history  of  loading  at  a particular  point,  since  for  almost 
all  of  the  elements  the  element  encompassing  that  point  will  be  the 
same  element  for  the  different  meshes  for  different  points  in  time. 


r 


3 


Figure  1.  Typical  Plug  Nose  Tip  Blocked  Into  Rectangular  Mesh. 


The  generation  of  the  mesh  proceeds  as  follows:  Given  the  complete 


outside  boundary  contour  definition  via  a series  of  points,  the  set 
of  reqular  rectangular  elements  encompassing  this  shape  may  be  determined 
by  simply  checking  the  intersection  that  each  R = constant  generator 
and  Z = constant  generator  makes  with  the  outside  boundary  contour. 
Knowing  the  location  of  the  intersections  along  these  generators, 
then  the  immediate  next  inter-element  location  is  chosen  as  the 
generator's  termination  point.  This  wilL  ensure  that  these  R = 
constant  and  Z = constant  generators,  which  are  simply  the 
lines  describing  the  boundaries  between  the  basic  rectangular 
elements,  describe  the  correct  distribution  of  rectangular  elements 
to  comprise  the  final  mesh.  Using  this  technique,  one  will  notice  in 
figures  1 and  2 that  this  initial  selection  of  rectangular  elements 
both  protrudes  out  from  and  into  the  area  defined  by  the  given  outside 
boundary  contour.  Hence,  the  next  task  is  to  adjust  those  rectangular 
elements  lying  next  to  or  on  this  outside  boundary  contour  so 
that  the  final  mesh  will  be  composed  of  reasonably  shaped  elements 
along  this  boundary.  However,  this  initial  sizing  using  rectangular 
elements  fixes  the  nodal  point  I-J  grid  for  the  mesh,  which  is  probably 
the  most  difficult  task  when  generating  a mesh  automatically  as 
described  here.  Further  information  about  the  relationship  between 
the  mesh  nodal  point  I-J  grid  and  its  R-Z  coordinate  grid  is  given 
in  Section  1 of  Reference  2 and  therefore  will  not  be  discussed  here. 

The  terminating  points  of  the  rectangular  element  arrays  comprising 
the  nose  tip's  mesh  as  shown  in  Figures  1 and  2 are  the  "boundary 
points"  of  these  meshes,  and  may  be  classified  as  one  of  three  different 
types  : type  T = triangular  point  type  C = corner  point  and 
type  E = end  point  (a  fourth  type  A = apex  point  is  a special  type 
and  applies  to  only  one  point  in  the  mesh) . These  different  types 
of  points  are  graphically  illustrated  in  Figures  3 and  4 , where  the 


symbols  T,  C and  E are  shown  next  to  the  type  of  boundary  point  which 
they  describe. 


A triangular  (type  T)  boundary  point  is  one  in  which  both  of  the  two 


sides  of  the  basic  rectangular  element  connecting  to  it,  intersect 


the  boundary  contour,  somewhere  along  these  two  connecting  sides. 


Consequently,  this  triangular  point  and  these  two  intersections 


comprise  a triangle;  hence,  the  name  triangular  point.  These  triangles 


are  shown  shaded  in  Figures  3 and  4 , and  as  shown,  can  exist  either 


outside  or  inside  of  the  boundary  contour. 


A corner  (type  C)  boundary  point  is  one  in  which  one  of  the  two 


sides  of  the  basic  rectangular  element  connecting  to  it  intersects 


the  boundary  contour,  and  the  other  connecting  side  doesn't.  Typical 


(C)  points  are  shown  in  Figure  3 . 


An  end  (type  E)  boundary  point  is  one  which  does  not  fall  into 


the  two  previous  types  of  classifications,  but  simply  represents  the 


"end"  of  either  a R = constant  or  Z = constant  generator.  Hence  the 


name  end  point.  Typical  (E)  points  are  shown  in  Figures  3 and  4 . 


The  readjustment  of  triangle  points  is  performed  as  follows: 


The  two  intersections  of  the  connecting  sides  define  a straight  line 


segment  comprising  part  of  the  structural  mesh  outside  boundary, 


shown  by  a dashed  line  for  a typical  triangular  point  in  Figure  5a. 


A perpendicular  is  formed  from  the  triangular  point  (T)  to  this  line 


segment,  and  the  intersection  of  this  perpendicular  and  the  (dashed) 


line  segment  comprises  the  final  location  of  the  triangular  point  on 


the  boundary,  as  shown  in  Figure  5a.  Hence,  the  triangular  point 


is  simply  translated  to  the  boundary  contour  via  this  "perpendicular' 


technique . 
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Figure  3.  Initial  Sizing  of  Nose  Section  of  Plug  or  Shell  Nose  Tip. 


Figure  4.  Initial  Sizing  of  Midsection  of  a Plug  Nose  Tip 


Final  Location  of 
Triangular  Point 
or  Boundary 


(a) . Triangular  Point. 


Figure  5(a)  and  (b) . Typical  Triangular  and  Corner  Points. 


The  readjustment  of  corner  points  is  performed  as  follows: 

The  one  connecting  side  of  the  basic  rectangular  element  which  intersects 
the  boundary  is  taken,  along  with  the  element's  parallel  side  (which 
also  intersects  the  boundary),  to  again  yield  two  intersections  which 
define  a straight  line  segment  comprising  part  of  the  structural 
mesh  outside  boundary.  This  is  shown  by  the  dashed  line  in  Figure  5. 
These  two  intersections  are  simply  averaged,  yielding  the  midpoint 
of  this  boundary  line  segment,  thus  comprising  the  final  location  of 
the  corner  point  on  the  boundary.  Hence,  a corner  point  is  simply 
translated  to  the  boundary  via  this  "averaging"  technique. 

The  readjustment  of  an  end  point  is  simply  performed  by  translating 
the  point  to  the  boundary  location  which  is  the  intersection  of  the 
boundary  and  the  R = constant  or  Z - constant  generator  which  the  end 
point  describes.  This  type  of  simple  translation  is  shown  for  end  (E) 
points  in  Figures  3 and  4 . 

The  fourth  special  type  of  point,  called  an  apex  point,  exists 
for  plug  type  of  nose  tips.  This  is  the  point  which  is  the  intersection 
of  the  outside  front  surface  and  the  ramp  portion  of  the  backside 
surface,  and  is  indicated  as  such  in  Figure  1.  This  point,  labeled 
(A)  in  Figure  4 is  simply  translated  to  the  actual  location  of  the 
apex  defined  by  the  input  data.  This  translation  is  shown  in  Figure  4 

After  all  of  the  triangular,  corner  and  end  points  have  been 
readjusted  to  lie  on  the  boundary  of  the  nose  tip,  the  original 
basic  rectangular  elements  connecting  to  these  points  become  general 
quadrilateral  elements  along  the  boundary,  as  indicated  by  dashed 
lines  in  Figures  3 and  4 . However,  some  of  these  quadrilateral 
elements  are  ill-shaped,  as  shown  by  the  insert  in  Figure  3 , i.e., 
two  of  the  four  sides  of  a general  quadrilateral  element  almost  form 
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a straight  line  (shown  dashed  in  the  insert  figure)  indicating  that 
this  particular  element  is  better  represented  by  a triangular  element 
than  a quadrilateral.  Hence,  all  of  the  boundary  quadrilateral  elements 
are  checked  for  this  type  of  ill-conditioning,  and  if  bad  enough  they 
are  changed  into  compatible  triangular  elements  by  simply  eliminating 
the  boundary  midpoint  (see  insert  in  Figure  3) . 


Thus,  the  structural  finite  element  mesh  of  a plug  or  shell 
nose  tip  is  automatically  generated  from  knowing  only  the  boundary 
contour  of  the  nose  tip.  A typical  mesh  of  a plug  and  shell  nose  tip, 
generated  by  PGMESH,  is  shown  in  Figures  13  and  15  respectively, 
where  the  small  "box''  symbols  plotted  along  the  boundaries  represent 
the  original  set  of  points  received  from  the  PAGAN  code  input  data, 
defining  these  boundaries. 


SECTION  III 


PGPRES  AUTOMATIC  NOSE  TIP  PRESSURE  FORCE  TRANSLATOR 

The  set  of  subroutines  contained  in  the  segment  PGPRES  are  almost 
identical  to  the  DOASIS  preprocessor  program  PRSINT,  described  in 
Section  4 of  Reference  2.  The  boundary  contour  for  both  the  input 
PAGAN  mesh  and  the  generated  structural  finite  element  mesh  are  already 
known  and  reside  in  core  when  PGPRES  is  called,  and  therefore  need 
not  be  input  as  described  in  the  description  of  PRSINT.  The  actual 
translation  technique  used  in  PGPRES  is  the  "normalized  arc  length 
translation"  technique  described  in  Section  4.3  of  Reference  2. 

The  backside  reacting  surface  is  described  exactly  the  same  way  as 
in  PRSINT,  where  the  variable  ISHPLG  defines  the  type  of  reacting 
surface  (see  Figure  20,  Reference  2) . 

Since  the  technique  and  theory  are  completely  described  for  the 
preprocessor  program  PRSINT  and,  therefore,  PGPRES,  in  Section  4 of 
Reference  2,  they  will  not  be  covered  here.  The  only  difference  between 
PRSINT  and  PGPRES  is  that  in  addition  to  translating  the  surface  pressures, 
PGPRES  also  simultaneously  translates  the  surface  temperatures,  thereby 
eliminating  a duplication  in  coding  of  the  surface  temperature  translation 
normally  included  in  the  segment  PGTEMP.  Since  the  surface  translation 
technique  is  the  same  for  pressures  in  PRSINT  and  temperatures  in  TEMINT, 
both  of  these  surface  translations  were  incorporated  in  PGPRES  instead 
of  including  them  seperately  in  PGPRES  and  PGTEMP. 
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SECTION  IV 


PGTEMP  AUTOMATIC  NOSE  TIP  TEMPERATURE  FIELD  TRANSLATOR 

The  set  of  subroutines  contained  in  the  segment  PGTEMP  are  very 
similar  to  the  DOASIS  pre-processor  program  TEMINT,  described  in 
Section  3 of  Reference  2,  with  one  major  difference.  That  difference 
is  locating  in  the  PAGAN  mesh  the  "element"  that  contains  the  point 
whose  temperature  is  being  translated.  Again,  no  need  exists  to  input 
either  the  PAGAN  mesh  or  the  structural  mesh  since  they  are  both  resident 
in  core  when  PGTEMP  is  called. 

The  input  PAGAN  mesh  contains  basically  two  types  of  elements: 

(1)  a rectangular  element,  basically  contained  in  the  interior  of  the 

mesh,  and  (2)  a rectangular  element  in  which  one  or  more  corners  are 

cut  off  by  the  boundary  contour,  i.e.,  transition  elements  lying  along 

the  boundary  of  the  PAGAN  mesh.  To  help  illustrate  this  point,  a 

typical  PAGAN  plug  nose  tip  mesh  is  shown  in  Figure  6.  One  can  observe 

that  the  rectangular  elements  lying  along  the  boundary  have  from  one 

to  three  corners  missing.  This  point  is  further  illustrated  in  Figure  7, 

which  shows  a rectangular  element  that  has  one,  two,  and  three  corners  missing. 

When  two  and  three  corners  of  a boundary  rectangular  element  are  missing, 
this  rectangular  element  simply  becomes  a general  quadrilateral  or  a 
triangular  element  respectively,  as  indicated  in  Figure  7.  However, 
when  only  one  corner  is  missing,  then  the  shape  of  the  rectangular 
element  is  a quintic,  that  is,  a polygon  having  five  vertices.  Inter- 
polation within  such  a polygon  could  be  performed  via  a Schwarz-Chr istof fel 
technique,  but  it  proves  easier  to  simply  subdivide  the  quintic  into  a 
triangle  and  quadrilateral  as  indicated  in  Figure  7.  This  subdivision  is 
performed  in  only  two  of  the  five  possible  ways.  In  both  ways  the  boundary 
line  segment  comprises  one  side  of  the  triangle  as  shown  by  the  top  two 
diagrams  in  Figure  7.  The  aspect  ratio  of  the  triangle  is  checked 
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for  these  two  possibilities,  and  the  one  with  the  best  conditioning 
is  chosen  as  the  subdivision  to  use  for  interpolation  purposes.  Hence, 
all  of  the  PAGAN  mesh  boundary  "transition"  elements  may  be  represented 
by  general  triangular  and  quadrilateral  elements. 

The  temperature  translation  for  the  interior  nodal  points  is 
performed  as  follows  (remember  that  the  boundary  nodal  point  temperatures 
were  interpolated  by  PGPRES) : 

1.  All  of  the  interior  rectangular  elements  are  checked  to  see 
if  the  nodal  point  lies  within  them,  and  if  one  is  found, 
it  is  used  for  the  interpolation. 

2.  If  none  of  the  interior  rectangular  elements  contains 
the  nodal  point,  it  must  lie  within  one  of  the  boundary 
transition  elements.  They  are  all  checked  and  the  one 
containing  the  nodal  point  (either  triangular  or  quadrilateral) 
is  used  for  the  interpolation. 

The  technique  and  theory  for  translating  temperatures  within  a 
general  triangular  and  quadrilateral  element  is  completely  described 
for  the  preprocessor  program  TEMINT  in  Section  3 of  Reference  2,  and 
therefore  will  not  be  reiterated  here.  This  applies  to  all  of  the 
transition  element  interpolations. 

The  technique  for  translating  within  a rectangular  element  is 
exactly  the  same  as  for  a general  quadrilateral,  described  in  Section 
3.3  of  Reference  2.  The  only  difference  being  the  determination  of 
the  non-dimensional  coordinates  £*  and  g*.  They  are  simply  calculated 
as  follows  (see  Figure  8). 

= 2 (f)  -1  g*  = 2 -1 


17 


Knowing  the  location  of  the  nodal  point  (£*,n*)  within  the  rectangular 
element,  the  temperature  is  interpolated  using  Equation  (24)  of 
Reference  2 via 

T(r,z)  = I (l-C)(l-n)T1  + <l+£)(l-n)T2 

+ (1+0  (l+n)T3  + ( l— £ ) (i+n)T4 
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SECTION  V 


PGP LOT  AUTOMATIC  NOSE  TIP  MESH  AND  TEMPERATURE  CONTOUR  PLOTTER 

The  set  of  subroutines  contained  in  the  segment  PGPLOT  are  almost 
identical  to  the  DOASIS  post  processor  program  CONTUR,  described  in 
Section  6 of  Reference  2.  The  plotting  software  used  is  identical  to 
that  described  in  Section  1.1  of  Reference  2,  and  is  reproduced  for 
convenience  here  in  Appendix  B.  The  mesh  plotter  and  temperature 
contour  plotter  are  the  same  as  that  contained  in  CONTUR,  with  minor 
differences  in  connection  with  the  title  actually  plotted. 

Since  the  technique  and  theory  are  completely  described  for  the 
post  processor  program  CONTUR,  and  therefore  PGPLOT,  m Section  6 of 
Reference  2,  they  will  not  be  covered  here. 

1 
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DESCRIPTION  OF  POESSY  COMPUTER  CODE 

A complete  listing  of  the  source  card  images  is  not  presented 
here  due  to  its  extensive  length.  Instead,  a block  outline  of  the 
different  subroutines  and  functions  is  shown  in  Figure  9.  The  subroutines 
PGMESH , PGPRES , PGTEMP  and  PGPLOT  shown  in  "dark"  outlined  boxes, 
call  other  subroutines  to  perform  their  functions  and  therefore  are 
shown  by  separate  block  outline  diagrams  in  Figure  9.  The  Plot  Software 
for  PGPLOT  is  shown  in  a "dashed"  outlined  box,  and  is  explained  in 
Appendix  B of  this  report.  There  are  three  tape  units  (which  may  be 
tapes,  disc  files,  drum  files,  or  any  other  peripheral  units  available) 
shown  in  Figure  9.  TAPE  2 is  the  PAGAN  code  data  input  tape.  TAPE  3 
is  the  save  tape,  used  to  store  (in  binary)  the  resultant  nose  tip 
finite  element  meshes  for  later  use  by  DOASIS.  TAPE  9 is  the  plot  tape. 

Table  1 is  a brief  summary  of  each  subroutine  and  the  function 
that  it  performs.  This  will  be  done  in  alphabetical  order  or  the 
subroutine  name. 


TABLE  1 


APEXEL 


AREA 


ASPECT 


BCELNO 


CONNEC 


CONTOR 


CONVAL 


SUBROUTINES  AND  FUNCTIONS 

Checks  the  "apex"  element  of  either  a shell  or  pluq  nose 
tip  to  see  if  a given  point  (P,Q)  lies  within  it. 

Calculates  the  area  of  a traiagle  given  the  coordinates 
of  it's  three  vertices. 

Checks  the  aspect  ratio  of  a triangle  and  determines 
the  ratio  of  the  height  to  the  length  of  the  longest 
side  of  the  triangle. 

Determines  the  boundary  contour  element  number  (M) , 
given  two  of  the  element's  nodal  points  lying  on  the 
boundary  of  the  finite  element  mesh. 

Generates  the  element  conductivity  information  by 
establishing  the  four  nodal  point  numbers  which  define 
the  element  in  the  finite  element  mesh. 

Generates  a isotherm  contour  plot  on  the  resultant 
translated  temperature  field  of  a nose  tip,  that  is,  it 
draws  lines  of  constant  temperature:  Either  automatic 

or  user  supplied  scaling  may  be  used  to  make  this 
temperature  contour  plot.  It  also  annotates  the  contour 
plot  with  symbols  along  the  boundary  contour  indicating 
different  values  of  constant  temperature. 

Interrogates  the  temperatures  being  plotted,  to  determine 
the  maximum  and  minimum  values,  and  then  automatically 
chooses  plot  values  of  the  temperatures  such  that  a 
minimum  number  of  values  exist  in  the  final  determination. 


TABLE  1.  Continued. 


CORNER 


DASHLN 


DAT  IN 


DOSAVE 


DRAW 


Eliminates  a "corner  point"  which  exists  on  the  boundary 

1- J  grid  to  help  smooth  out  the  final  boundary  surface 
contour  in  a manner  similar  to  the  subroutine  TRIANG. 

Draws  a dashed  line  similar  to  the  plotting  software 
routine  LINE. 

Inputs  the  nose  tip  mesh  information  from  a PAGAN  code 
save  tape,  either  via  reading  from  cards  or  reading  from 
the  PAGAN  code  save  tape.  It  redistributes  this  input 
information  into  a form  which  may  be  used  by  the  POESSY 
code,  that  is,  it  redistributes  the  coordinate  and  I-J 
grid  information  into  the  POESSY  internal  format. 

It  redistributes  the  interior  coordinate  points  into 

2- dimensional  arrays  which  are  used  for  mesh  generation 
and  temperature  translation  purposes.  It  checks  the 
list  of  surface  nodal  point  numbers  to  see  if  any  surface 
nodal  point  is  duplicated  or  belongs  to  the  interior 
arrays,  and  eliminates  any  duplications. 

Writes  the  final  structural  finite  element  mesh  onto  the 
DOASIS  save  tape  in  a format  compatible  for  inputting  into 
the  DOASIS  finite  element  program. 

Draws  contour  plots,  given  the  coordinate  variables  and 
the  function  variable  in  2-dimensional  arrays  and  the 
appropriate  maximum/minimum  limits  of  these  2-dimensional 
arrays.  Additional  features  included  the  ability  to 
exclude  certain  function  variable  values  from  the  final 
contour  plot. 
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TABLE  1.  Continued. 


ELEMCK 


ELEMNO 


ELIMTP 


ELLABL 


ELMS I Z 


Checks  to  see  if  a given  point  (P,Q)  lies  within  or  on 
a quadrilateral  element  by  subdividing  the  quadrilateral 
into  four  triangles,  formed  by  the  point  (P,Q)  and  one 
of  the  quadrilateral’s  sides.  It  then  checks  to  see  if 
the  area  of  these  four  triangles  is  positive,  and  if  so, 
establishes  that  the  point  lies  within  the  given  quadrilateral 
element . 

Determines  the  element  numbers  of  elements  which  connect 
to  a given  nodal  point  (N) . 

Checks  the  degree  of  singularity  which  might  exist 
at  a corner  of  a surface  quadrilateral  element,  and  changes 
this  quadrilateral  element  into  a compatible  triangular 
element  if  the  singularity  exceeds  a predetermined 
(user  supplied)  criteria. 

Annotates  the  boundary  of  the  finite  element  mesh  with 
the  list  of  boundary  contour  element  numbers,  placed 
appropriately  next  to  the  elements  on  the  boundary. 

Determines  the  basic  rectrangular  element  size  (DI,DJ) 
based  upon  the  total  area  of  the  nose  tip  being  composed 
of  approximately  (NELEM)  elements.  If  either  or  both 
of  the  basic  spacing  data  DI  and/or  DJ  is  given,  then 
these  values  are  used  when  determining  the  basic 
rectrangular  element  size.  If  neither  are  given,  then 
these  values  are  choosey  based  upon  the  nose  tip  being 
composed  of  NELEM  elements. 
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TABLE  1.  Continued. 

FETCHA  --  Extracts  an  alpha-numeric  character  from  a computer  word, 

given  the  position  in  the  word. 

FILBAK  --  Adds  to  the  list  of  backside  boundary  surface  nodal 

point  numbers,  those  points  that  were  given  in  the 
interior  definition  of  the  mesh  and  were  not  originally 
included  in  the  complete  backside  boundary  list. 

FIND  --  Locates  the  point  in  the  PAGAN  mesh  where  a given  nodal 

point  from  the  structural  finite  element  mesh  lies,  and 
then  translates  the  temperature  of  this  given  structural 
nodal  point.  It  first  checks  all  of  the  basic  PAGAN 
rectangular  elements  to  see  if  the  point  lies  within 
one  of  these  elements  and,  if  so,  it  uses  this  basic 
rectangular  PAGAN  element  to  translate  the  temperature. 

If  the  point  does  not  lie  within  one  of  these  basic 
rectangular  elements,  then  all  of  the  transition  elements 
existing  along  the  outside  and  backside  surfaces  are 
checked  and  the  subroutines  TRI  or  QUAD  are  called  to 
translate  the  given  temperature. 

FTXMM  --  Adjusts  the  nodal  point  I-J  grid  maximum/minimum  information 

to  reflect  the  fact  that  a boundary  contour  nodal  point 
has  been  eliminated  from  the  mesh  definition. 

FIXPNT  — Finds  the  intersection  of  two  lines,  each  of  which  is 

defined  by  two  nodal  point  numbers. 

FIXSPC  --  Changes  the  units  and  reference  datum  of  the  basic  spacing 

data  received  from  the  PAGAN  save  tape,  that  is,  the  basic 
spacing  variables  RSPACE  and  ZSPACE. 
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TABLE  1.  Continued. 


FIXSUR 


FORCE 


FPNUM 


FRAMIT 


GETIR 


Checks  the  surface  nodal  point  number  lists  to  see  if 
any  have  been  flagged  for  elimination  and  eliminates 
these  points  from  the  final  surface  nodal  point  list. 

It  reorders  this  list  of  nodal  point  numbers  to  reflect 
the  fact  that  certain  points  have  been  eliminated  and 
updates  the  nodal  point  I-J  grid  maximum/minimum  arrays 
to  reflect  that  certain  boundary  contour  nodal  points 
have  been  eliminated  from  the  resultant  structural 
finite  element  mesh  definition. 

Calculates  the  total  axial  and  radial  forces  existing 
on  a boundary  surface  line  segment  which  is  subjected 
to  normal  and  tangential  pressures.  It  assumes  that 
these  pressure  forces  are  constant  along  this  line  segment. 

Encodes  a floating  point  number  from  internal  format  into 
equivalent  alpha-numeric  (F)  format  for  plot  annotation. 

Automatically  determines  all  of  the  necessary  sizing  data 
to  make  the  resulting  plots  fit  into  the  specified  paper 
size,  that  is,  the  paper  length  and  the  paper  width, 
given  the  parameters  controlling  the  plot,  such  as,  the 
title  information,  the  margin  specifications,  the  rotation 
and  scaling  specifications  and  the  maximum  and  minimum 
coordinate  information. 

Determines  the  basic  spacing  index  (IR)  corresponding 
to  a given  radial  coordinate  value  belonging  to  the  list 
of  radial  coordinate  spacing  values  (RSPACE) . 
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TABLE  1.  Continued. 


GETJZ 


IJBCNO 


IJLIM 


IJSIZE 


INNUM 


IRJZCK 


KLOCK 


Determines  the  basic  spacing  index  (IZ)  corresponding 
to  a given  axial  coordinate  value  belonging  to  the  list 
of  axial  coordinate  spacing  values  (ZSPACE) . 

Generates  the  element  I-J  grid  maximum/minimum  information 
from  the  nodal  point  I-J  grid  maximum/minimum  information. 

Determines  the  closest  basic  spacing  data  index  I or  J, 
given  the  value  of  a coordinate  and  the  basic  spacing 
dimension  of  that  coordinate,  i.e.,  given  R and  DI  or  Z 
and  DJ. 

Initially  determines  the  I-J  grid  information  associated 
with  the  intersections  of  the  fixed  DI-DJ  grid  of  the 
structural  mesh  and  the  outside  and  backside  surfaces 
defined  by  the  PAGAN  input  data.  This  determines  the 
initial  I-J  grid  maximum/minimum  information  associated 
with  the  structural  finite  element  mesh. 

Encodes  an  integer  number  from  internal  format  into 
equivalent  alpha-numeric  (I)  format  for  plot  annotation. 

Checks  to  see  if  a point  defined  by  the  indical  values 
(IR,JZ)  lies  outside  of  the  PAGAN  mesh  IG-JG  grid. 

Calls  a system  subroutine  to  determine  the  current 
elapsed  central  processor  time. 
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TABLE  1.  Continued. 


t 


KPLACE  — Determines  the  location  and  word  of  a character  for 

encoding  integer  and  floating  point  numbers,  as  required 
by  the  subroutines  INNUM  and  FPNUM. 

LIMITS  — Determines  the  maximum/minimum  values  of  the  two  coordinate 

arrays  X,Y. 

MAKTRI  — Changes  a quadrilateral  element  into  a compatible 

triangular  element  by  renumbering  the  element's  nodal 
point  numbers  to  reflect  the  fact  that  one  of  them  has 
been  eliminated. 

MINMAX  --  Determines  the  maximum  and  minimum  I-J  grid  information 

for  the  structural  finite  element  mesh. 

MOUTCK  --  Checks  to  see  if  a element's  (IC-JC)  indices  lie  outside 

of  the  element  I-J  grid,  thereby  indicating  that  this 
element  does  not  exist  in  the  structural  element  mesh. 

MXMN  — Establishes  the  PAGAN  mesh  I-J  grid  maximum/minimum 

index  information. 

NMIJNP  — Performs  the  actual  function  of  checking  the  boundary 

surface  elements  to  see  if  a given  point  (P,Q)  lies  inside 
one  of  them. 

NODE  — Determines  the  nodal  point  number  (N) , given  its  indices 

(I,J)  on  the  nodal  point  I-J  grid. 


31 


TABLE  1 . Continued . 


NODE  I J 


NOELEM 


NORMAL 


NOUTCK 


NPLABL 


NPLOC 


NPPRFS 


Determines  a nodal  point's  indices  (I,J)  on  the  nodal 
point  I-J  grid,  given  its  nodal  point  number  (N). 

Returns  the  element  number  (M)  given  its  indices 
(IC,JC)  on  the  element  I-J  grid. 

Calculates  various  geometrical  quantities  pertaining  to 
the  structural  surfaces,  e.g.,  surface  line  segment 
lengths,  angular  orientations  and  outward  normal  directions 
for  both  line  segments  and  surface  nodal  points. 

Checks  to  see  if  a nodal  point's  (I,J)  indices  lie 
outside  of  the  nodal  point  I-J  grid,  thereby  indicating 
that  this  nodal  point  does  not  exist  in  the  structural 
mesh . 

Annotates  the  boundary  of  the  finite  element  mesh  with 
the  list  of  boundary  contour  nodal  point  numbers,  placed 
appropriately  next  to  the  nodal  points  on  the  boundary. 

Determines  the  plot  coordinate  location  of  the  number 
annotation,  required  by  the  subroutines  NPLABL  and 
ELLABL. 

Generates  the  backside  surface  definition,  that  is, 
the  nodal  point  numbers  and  the  associated  coordinate 
values  of  a plug  or  shell  nose  tip.  The  particular  type 
of  backside  surface  definition  is  determined  by  the 
user  via  an  input  option. 


TABLE  1. 


Continued . 


NSBLIM 


NSURNP 


PGMESH 


PGPLOT 


Searches  the  boundary  contour  nodal  point  lists  to 
determine  the  location  and  number  of  outside  and  backside 
surface  nodal  point  numbers. 

Checks  to  see  if  a given  nodal  point  number  (N)  lies 
on  the  boundary  surface,  that  is,  whether  it  exists  in 
the  boundary  surface  nodal  point  number  lists. 

Main  subroutine  which  controls  the  automatic  generation 
of  nose  tip  mesh.  It  calls  subroutines  to  input  the 
PAGAN  data  from  a save  tape  and  to  automatically  size 
the  nose  tip  for  a given  element  size.  It  sets  up  the 
nodal  point  I-J  grid  for  the  nose  tip,  and  calculates 
the  maximum/minimum  I-J  grid  information  associated  with 
the  nose  tip.  It  fills  in  the  interior  R and  Z coordinates, 
for  the  nose  tip  and  determines  the  element  I-J  grid  for 
the  nose  tip.  It  generates  the  outside  surface  boundary 
contour  nodal  point  lists.  It  checks  for  bad  quadrilateral 
elements  which  might  exist  along  the  outside  boundary  of 
the  nose  tip,  and  if  needed,  smooths  out  these  bad 
quadr ilaterials  by  changing  them  into  triangular  elements. 
Finally,  it  prints  out  the  results  of  the  mesh  generation. 

Main  subroutine  which  controls  the  plotting  of  the  final 
structural  finite  element  mesh  generated  by  PGMESH . 

It  sizes  the  plot  to  fit  onto  a specified  paper  size, 
and  also  annotates  the  plot  with  either  boundary 
contour  nodal  point  or  element  numbers.  It  also  plots 
the  temperature  field  which  was  translated  by  PGTEMP  by 
calling  CONTOR. 
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TABLE  1. 


Continued . 


PGPRES 


PGTEMP 


PLTFN 


PLTST 


POESSY 


PSCALE 


Main  subroutine  which  translates  the  nose  tip  pressure 
distribution  given  by  the  PAGAN  input  data  to  the  generated 
structural  finite  element  mesh.  It  translates  the  boundary 
pressures  and  temperatures  for  the  outside  surface  boundary 
and  calculates  a reacting  pressure  distribution  for  the 
backside  of  the  nose  tip  such  that  the  resultant  pressures 
acting  on  the  nose  tip  are  in  static  equilibrium. 

Main  subroutine  which  translates  the  PAGAN  mesh  temperature 
field  to  the  resultant  structural  finite  element  mesh. 

This  translation  process  is  performed  by  determining  the 
temperature  of  the  nodal  points,  and  then  averaging  these 
nodal  point  values  to  determine  the  element  temperatures. 

Terminates  the  plotting  sequence  for  a given  set  of  plots. 

Initializes  the  plotting  routines. 

Main  driving  program  which  calls  the  four  sub-programs 
PGMESH , PGPRES,  PGTEMP,  PGPLOT,  to  perform  the  mesh 
generation,  pressure  translation,  temperature  translation, 
and  plot  the  resulting  mesh  and/or  temperature  contours, 
respectively. 

Automatically  scales  the  axes  of  the  plot,  given  the  axis 
length  of  one  coordinate  direction  and  the  maximum/minimum 
information  of  both  coordinate  directions. 
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QDRANT  — Determines  the  quadrant  in  which  a line  defined  by  the  two 

incremental  lengths  DR  and  DZ,  lies. 

QUAD  --  Interpolates  the  temperature  of  a point  lying  within  a 

general  triangular  element. 

QUADSZ  — Determines  the  aspect  ratio  associated  with  a quadraliteral 

finite  element. 

RCHECK  — Checks  to  see  if  a given  radial  coordinate  R is  equal 

to  one  of  the  basic  spacing  data  points  RSPACE. 

RESULT  --  Prints  (and  punches  in  DOASIS  compatible  format)  the 

resultant  structural  finite  element  mesh. 

ROUND  — Rounds  a given  floating  point  value  to  a specified  number 

of  digits. 

RZPRIN  — Prints  out  the  structural  finite  element  mesh  grid  information 

that  is,  the  radial  and  axial  coordinates  in  2-dimensional 
array  form,  the  nodal  point  and  element  I-J  grid  maximum/ 
minimum  information,  the  list  of  boundary  contour  nodal 
point  numbers  and  the  list  of  outside  and  backside  surface 
nodal  point  numbers  along  with  their  corresponding  coordinate 
values . 

SBIJMN  — Determines  the  I and  J indices  associated  with  the  basic 

spacing  data,  RSPACE,  ZSPACE  and  the  limits  of  the  R and 
Z coordinates. 


TABLE  1.  Continued. 


k 


SLOPE  --  Determines  the  constants  A,  B,  C and  D of  the  equations 

defining  a line,  that  is,  Z = A*R  + B,  or  R = C*Z  + D, 

when  given  the  incremental  values  DR,  DZ  and  RZ . 

SMOOTH  --  Checks  and  smooths  out  the  irregular  boundary  surface 

coordinates  produced  by  the  original  sizing  of  the  nose 
tip.  Points  which  lie  outside  of  the  original  surface  are 
moved  to  the  boundary  surface  in  such  a manner  that  the 
resulting  quadrilateral  and  triangular  elements  lying  along 
the  surface  form  a smooth  and  even  array  of  elements. 

Points  that  lie  inside  of  the  original  boundary  surface 
are  similarily  moved  to  the  boundary  surface  locations  to 
create  meaningifully  shaped  elements  comprising  the  final 
structural  finite  element  mesh. 

SORTR  — Sorts  the  R and  Z coordinate  lists  and  (T)  temperature  list, 

according  to  decreasing  values  of  R. 

STOREA  --  Stores  an  alpha-numeric  character  into  a computer  word, 

given  the  position  in  the  word. 


SURBAK 


SURFEL 


Extracts  from  the  boundary  contour  nodal  point  lists  the 
list  of  outside  and  backside  nodal  point  numbers  and  coor- 
dinates associated  with  the  outside  and  backside  pressure 
surface  definitions,  and  loads  the  appropriate  arrays  with 
these  numbers  and  coordinates. 


Checks  all  of  the  boundary  surface  elements  of  the  PAGAN 
mesh  to  see  if  a given  point  (P,Q)  lies  inside  or  on 
one  of  them.  If  so,  it  calls  the  appropriate  subroutine 
to  find  the  temperature  of  this  given  point  (P,Q). 
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TABLE  1.  Concluded. 


SURPIJ 


SYMLAB 


TITPLT 


TRI 


TRIANG 


TRISIZ 


ZCHECK 


Performs  three  distinct  functions;  (1)  it  checks  R and  Z 
coordinate  arrays  and  eliminates  any  points  from  these 
strings  which  are  duplicates;  (2)  it  checks  the  I and  J 
indices  associated  with  the  surface  nodal  points  to  see 
if  they  are  all  included  in  the  original  strings,  that  is, 
to  verify  that  all  of  the  intersections  of  the  interior 
mesh  and  the  outside  surface  boundary  are  present;  and 
(3)  it  assigns  the  I and  J indical  values  to  these  surface 
nodal  point  strings  which  are  used  later  in  the  solution 
sequence . 

Annotates  the  boundary  contour  with  temperature  value 

plot  annotations,  that  is,  symbols  placed  along  the  boundary 

which  represent  the  different  contour  temperature  values. 

Plots  title  information  on  the  plots,  relative  to  the 
origin  of  the  plot  axes  and  rotation  direction. 

Interpolates  the  temperatures  of  a point  lying  within  a 
general  triangular  element. 

Eliminates  a "triangular  point"  which  lies  either  inside 
or  outside  of  the  boundary  contour,  in  order  to  smooth 
out  the  resulting  finite  element  mesh  boundary  surface 
elements . 

Determines  the  aspect  ratio  associated  with  a triangular 
finite  element. 

Checks  to  see  if  a given  axial  coordinate  Z is  equal  to  the 
basic  spacing  data  points  ZSPACE. 
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SECTION  VII 


MODIFICATION  OF  PROGRAM  SIZE 

The  size  of  the  POESSY  computer  program  is  determined  by  three  factors. 
First,  the  coding  of  the  program  which  contains  all  locally  dimensioned  variables 
is  basically  fixed  by  the  computer  and  compiler  being  used.  Secondly,  the 
computer  being  used  also  determines  the  size  of  the  FORTRAN  library 
functions  and  system  routines  needed  to  run  POESSY.  Thirdly,  the  size 
of  the  common  blocks,  used  to  store  the  solution  variables,  is  determined 
by  the  user,  depending  upon  the  size  of  problems  he  wishes  to  solve. 

All  of  the  common  blocks,  and  relevant  data  statements  and  dimension 
and  equivalence  statements  are  shown  in  Figure  10,  along  with  the 
variables  ISIZ,  . . . , SSIZ.  In  order  to  change  the  size  of  the  code, 
one  simply  has  to  change  the  size  of  these  common  blocks,  dimension 
and  equivalence  statements  in  every  routine  in  whicn  they  are  contained. 

A cross-reference  table  to  SUBROUTINE  versus  COMMON  BLOCK  is  given 
in  Figure  11.  Therefore,  the  total  size  of  core  needed  for  all  the 
variables  contained  in  the  common  blocks  is  given  by 

LABELED  COMMON  = 1128  + 4*ISIZ  + 4*JSIZ  + NBSIZ  + 3*IGSIZ 

+ 3* JGSIZ  + 6*MELSIZ  + 4*NPSIZ  + 5*PSIZ 
+ 10*SSIZ  + 2*ISIZ* JSIZ  + 3*IGSIZ*JGSIZ 
where  the  meaning  of  ISIZ,  . . . , SSIZ  is  given  in  Figure  10. 

For  example,  if  the  sizing  is  that  shown  in  Figure  10,  then  the  sizes 
of  core  would  be 

LABELED  COMMON  = 22,478 
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The  size  of  core  needed  for  the  code  and  system  routines  is  given  here 
for  a Control  Data  Corporation  CDC  6600  (Extended  FORTRAN  Compiler) 
running  under  the  KRONOS  2.1  operating  system,  with  core  size  for 
common  blocks  repeated  from  above. 


CODE 

PLOT  CODE 

COMMON 

SYSTEM 

TOTAL 

CDC  6600 

23,262 

2,191 

22,478 

7,537 

33,053 

Similar  total  core  requirements  may  be  expected  from  other  computers. 
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Figure  11.  Cross  Reference  of  Subroutines  Versus 
Labeled  Common  Blocks  for  POESSY 
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Figure  11  • Concluded 
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I APE  5 A OUTPUT  DAT*  Tape  from  PTESsy  A InPuT  data  Tape  EuP  DUASIS 
TAPES  A T P*t  STAMjAHO  Input  EUt  E L P RtADING  CAPOS 
TAPES  = The  STANDARD  OUTPUT  EIlE  FuR  PKIM[N|, 

T A PE  7 A CAPU  PUNCH  PILE 

TAPES  £ OUTPUT  TAPE  E Op  PLOTTING  (uPUM  T YPf  PCOTTEPS) 


I'M.  EiJESST  CODE  ui.it  REDuIpES  one  CLnTkiil  C*pD  EUR  F Af  it  SETUP  DESIpEd. 
anu  Tali  ADI  IIIONAL  PLOT  SCALL1NG  CAPuS  IF  PLOTS  APE  hAuE . THESE  TmPEE 
C‘pOS  APt  DtSOIbfL  t>T  THE  EuLlOInG. 


* CUNlkUL  Cap u E UPMA T 1 9 1 5 , 3E I u . u ) 


COLUMNS  1 -S  MILASE  r CASE  I Tit  N T I E 1 1 A 1 ] ON  NUPEiER. 

b-10  ]ShP|_6  a n,  PLUG  ITPE  UE  NOSE  TIP,  «hEpE  ThE  REACTING 

PRESSURE  IS  ASSUMED  ACTING  on  Tht  BUM  EnO 
POP  T 1 C N UF  IhE  Shank,  IN  TnE  7 -D 1 PE C I I On . 

: 1,  SHELL  TYPE  OF  NOSE  TIP.  ahERE  THE  REACTING 
PRtSSUPfc  IS  ASSUMED  TU  ACT  ON  THE  INSIDE 
NOSE  PURIIUN  OF  The  SPtCL.  FROM  THE  AXIS 
OF  Wt  VUL  U T I On  to  a point  neap  IHt  SPhEPf* 
CONE  InTEpSECTIUn  PCInT  on  ThE  InSIDE, 

£ ?,  SHtLL  TYPE  Of  NOSE’  TIP.  "HERE  iMt  REACTING 
PRtSSUPE  IS  ASSUMED  TO  ACT  ON  IHt  BUTT  t ND 
OF  THE  CONE  PORTION  OF  ThE  SHELL  UP, 

11-1S  NSTEP  £ TIME  STEP  number  OF  Tht  PAkTICULAP  PuInT  In 
The  TPAJtClUPY  EuP  ahRm  TmE  me SH  SEIUP  IS 
DESIPFO  (This  IS  GIVEN  FROM  THt  pagan  SOLUTION) 
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PAGE  Nu 


i 


16-PO  I*MT  = PRIM  UUI  CLMROL  PARAMETER.  I Ht  RLtSSY  CODE 
MAS  SIX  (fe)  DIFFtRtNl  LEVELS  UE  PRINT  UUI 

available,  Ranging  fkum  IrKIT  = 0 TO  IwRIT  i 5, 

Ihtwt  MEANING  IS  BRIEFLY  AS  FOLLOWS... 

= U,  P R J f,  T II. Pul  CuNIkUL  iNFukMAlIUN  AND  ThE 
Summary  UE  [ME  UU1PUT  CUNIRuL  RESULTS, 

= 1,  ADDITIONALLY,  PRIM  OUT  THE  NODAL  PUInI 
COORDINATES,  TMt  ELEMENT  LONNECTIvITy 
INFORMA  I lUN,  TMt  PktSSURE  CARD  INFORMATION 
AND  THE  ELEMENT  T E MPE R A T UkE S , 

= 1?,  AoOITIOnALLY  PRIM  OUT  THE  1ST  LEVEL  OF 

ADDITIONAL  u I AGON  1 S T 1 C DATA  PERTAINING  TO 

INTERMEDIATE  calculations, 
i i,  ADDITIONALLY  prim  UUT  the  «>nd  level  of 

ADDITIONAL  DlAGUMSTIL  DATA  PERTAINING  to 
IMtRMtUlATE  CALCUl  ATIOnS. 

= R,  ADDITIONALLY  PRIM  OUT  TMt  3RD  LEVEL  uF 

additional  uugonistic  data  pertaining  to 

INTERMEDIATE  CALCULATIONS, 

= S,  AuDITiONALLY  PRINT  UUT  the  HlGHtST  LEVEL 
ADDITIONAL  DIa&UNISTIL  data  PERIAINING  to 
INTERMEDIATE  CALCULATIONS, 
ai-PS  1 PUNCH  r <1,  DU  NUT  PUNCH  UuT  CARDS. 

s 1,  PUNCH  UUI  The  mESm  NUOAL  POINT  CaRUS,  TMt 
ELEMENT  CONNECTIVITY  cards,  the  PRtSSDRE 
CARDS  AN D THE  1-J  N.P,  AND  ELEMtM  GRID 
CARDS, 

pfe-iU  IPLOT  = 0,  DU  NOT  PLOT  ANYTHING, 

= 1,  PLOT  THE  FINITE  ELEMENT  MESH, 

= i,  PLOT  THE  ISUTmERMS  (CUMOUR  PLOT), 
s i,  PLOT  B01H  THE  FINITE  ELE  Ht  N I mESH  And 
THE  I SO  I Ht  RMS  (CONTOUR  PLOT), 

J 1 -iS  IUSAVE  S 0,  SAVE  NOTHING  On  THE  UA1A  SAVE  TAPE  = TAPES. 

= 1,  save  the  complete  mesh  description  on  the 
save  TAPE  S TAPEi,  In  Ihe  DOASIS  CUDE 
INPUT  format, 

ib-RO  NELtM  s NUMBER  OF  ELEMENTS  Tu  DIVIDE:  NuSt  TIP  1M0 
WHEN  THF  BASIC  ELEMENT  SWt  IS  DETERMINED 
AUTOMATICALLY  (SEE  DiSIZ  AND  DJSIZ),  The 
FINAL  NUMRtR  OF  ELEMENTS  CANNU1  EXCEED  600, 
HENCE  A 1 YPICAL  VALUt  = SSO,  IF  NOTHING  IS 
SPECIFIED,  A Dtt  AUL  T VALUt  i R00  IS  USED. 

Rl-RS  1NPA&N  = 0,  INPUT  PAGAN  LODt  DATA  F RUM  TAPE?  = PAGAN 
CODE  S A V b TAPE  (NORMAL  OPTION), 
s 1,  INPUl  pagan  tout  DATA  FROM  CARDS  (TAPES), 
R6-5S  DISIZ  = BASIC  SPACING  DIMENSION  FOR  IHt  t Lt Mt N I SIZE 
IN  The  P-CuLRDlNAtt  DIRECTION, 

S6-6S  UJSJZ  = BASIC  SPACING  DIMENSION  FOR  THE  ELEMENT  SIZE 
IN  THE  Z-CUtJRD  I N A T t DIRECTION. 

*****  NOTICE  **»**,  ONE  MAY  SPECIFY  EIThfH 
DISIZ  AND/OR  DJSIZ  TO  PREDETERMINE  THE  BASIC 
ELEMENT  SIZING  INFORMATION,  IF  NtITHtR  IS 
SPECIFIED,  THtN  TmEy  ARt  DETERMINED  BY  THE 
PuESSY  CODE  HASEO  UPON  DIVIDING  TmE  NuSt  TIP 
JNTu  l(NELEM))  t L t ME  NT  S , AUTOMATICALLY. 
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PAGt  NO 


i 


00-7S  I R 1 E c I z UUAuRILAltRAL  ELEMtNI  RtjECIlUN  EaLIUr,  UStD 

TO  DElERMlNt  JF  IhusE  ELEMENTS  L'lNo  ALUM,  The 
BOUNDARY  lONTUUK  iiAvE  BAO  ASPtCT  RAlIUS,  AND 
SHOULD  HE  CHAuGtD  l N I 0 IhIanGuLAR  ELtMENls. 
(THE  DEEAULT  value  OF  IRIEL1  z -lib. 00) 

IF  TRIFLT  .01  . 0.0,  fnEN  1RIECT  RtPRtSENT& 

THE  MAXIMUM  ALLOWABLE  HtiOMl  10  BA  St  KAUu 
OF  A TRIANGLE  EURMtu  f HUM  THE  I wU  BAD  SIDtS  OF 
(JUAUH1LATEHAL  ELtHtM  . (TyRICAL  VALUE  = 0,1) 

IF  TRIELT  ,Ll.  0,0/  THt«-  l«IFLT  R t P Rfc  St  N T S 
THE  MAXIMUM  ALLOWABLE  fcNLLUStU  ANoLt  t A 1 S I I NO 
AT  T Ht  NODAL  F>UINT  CUNNtC  I I NO  Tu  T Ht  l»0  BAu 
SIDES  OF  IHt  ELEMtNI.  ( T VP  I L AL  vALUt  z lib.oO) 


* PLUT  SCALE  CUNTRUL  CAKU  FORM A T ( 8 lb/ bF 1 0 . 0 ) (UnLT  IE  IPlUI  .01,  0) 


THIS  CARD  CUNTROLS  THE  BoTATIUN  AND  SCALINO  OF  I HE  HLUI,  EIlHtR  FULLY 
AUTOMATIC  PROGRAM  LONTHUL  UR  USER  SuPPLitU  CUN T kUL  mat  Bt  UStD, 

COLUMNS  1-b  ISCALE  = 0/  T Ht  PROGRAM  will  AUTUMA I IC ALL ' SLALt  IHt 

CUUKDINAIt  DATA  10  FIT  THE  ASSUMED  P A Pt  M 
DIMENSIONS  (PAPtRL  AND  PAPtRw), 

= 1,  T Ht  UStK  »ILL  SUPPLY  THE  APPKUPlATt 

SCALING  INFORMATION  CRSTARI , ZS I A k T .OtLPO) , 
0- 1 0 IRUTAT  z -1,  PLOT  axes  WITH  Z - X-AXIS  UIKtLTION 
(PAPER  LENGTH  DIRECTION  z PAPERL) 
z 0 i MUTATE  I Ht  AXES  BASED  UPON  IHt  LONGER  OF 
ThE  TwO  AxES  GOING  IN  the  longer  of  Tht 
T wO  PAPER  DIMENSIONS  (LtNOTn  AND  wIuTm) 
z *1,  PLOT  AXES  h 1 Tm  R z x-AxlS  DIRECIlUN 
(PAPtR  LE  No  I H DIRECTION  z PAPtRL) 

11-lb  IMARON  z 0,  DU  NUT  Lt  A VE  A 1 INCH  MARGIN  HtTwttN  ThE 
BODY  OF  T Ht  STRUCT  uRt  AND  IHt  AXtS, 
z 1,  LEAVE  A 1 1 nlh  MARGIN  BETWEEN  THtM, 

10-20  IMXMN  r 0,  IGNORE  R S I AK T , l S T A R I , UELPu  (FuR  ISCALt  z 0) 
z 1,  USE  RSTART  AS  R-LOURD1 N A T t MINIMUM  value, 
z 2,  USE  ZSTARI  AS  Z-CUIIRO  IN  A Tt  MINIMUM  VALUt. 
z i,  USE  BuTH  RS I AR T , ZS 1 AH  I AS  MiNIMuMS. 

( R M I N AND/OR  ZM I N ARE  USED  IF  RSIAR1  AND/OR 
ZSIART  Awt  NUT  SPECIF ItO) 

2 1 -2b  InbTyP  z o,  l)U  NOT  AnnuTaTE  mesh  PLOT. 

z |,  anno  I A 1 E The  mesh  PluT  wITh  Ihe  BUUNDARY 

CON  I OUR  NUDAL  POINT  NUMBERS,  PLACEU  Nt X I TO 
iHtIR  LLCATIONS  OUT  SUlt  OF  ThE  MtSn, 

I 2,  ANN01AIE  IrE  MtSH  plot  wIIh  IHt  BOUNDARY 
ELt ME N 1 NUMBtRS,  PLALED  NEXT  To  IHtiR 
LOLAIIoNS  uulSIDt  HE  The  MtSH, 
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"tlLtk  PEStAkLk  InC.y  HI.  Vit*,  CALlf,  PAGE  nu.  h 


dO-iQ  1UHGLN  s Ur  Gtl  NUl  PLUI  Ikt  I NPU  I BIJUNGAkY  tuMlnjh 
GtEI'iIIlLN  I HUM  Ikt  HAGAN  LUGt, 
s 1.  PLUI  Tht  INPuI  BUUNgAHY  LUnIiJUh  Gt  PInIIIun 
F HUM  Ikt  Hagan  tHDt.  Ab  A gASheG  1-iNt.  »iln 
bVHrtULb  AT  Ikt  I - J Gkju  I N 1 E «5t C I 1 UNo , 
il-JS  1ANNUT  * 1.  fMt  VALUtS  liE  Ikt  CUMLUkb  PLUIltl-  »1LL 

bt  LlSlLG  UN  PLUl  y 

s Ur  Gu  Nt  I Llbl  Any  CUNlUuk  1NEU,  U I nt  PLuIy 
i -n,  IN  AGgIIIun  lu  lIoIINg  Ikt  vALUtS  UE  Ikt 
CUNIlukb  ON  Int  ktl  i.  A blttlAL  Symbll  | s 
USEG  TL  ANNuT  A I t LAtH  (N-lkJ  1 1 N T UGk  ALUNl, 
I kt  BUGNGAk  y tit  I Ht  BUG  V t ang  Ib  At  bU 
LISltG  Nt  * I It  Ikt  L IN  T UGk  v At  Ut  lArlt. 
5B-80  NCUNI  a HJnIHUH  NUMBtk  ut  ibUikfkkb  (ttMPtHAIUkt 

CIJMuUkS)  IU  PLUl.  Ikt  PUtSS*  tuGL  GtltkMJNtb 
I ht  FINAL  NGHot  H AnG  value  5 tit  Ikt  CUMUUkS 
Tl)  t>LUl  At)  I UHA  I It  ALL  Y . hENLEy  Ikt  FINAL  NIIMbtk 
Ut  LUNHukb  *ILL  Bt  .Gt*  1 u NtUNT  ANg  .Lt. 
d * N t UN  I . (UttAGtl  VALUE  IS  N t t)N  1 = lot. 


* HLU  I SLALt  LUNlkUL  C A k G . LUNlINUtl  l St  t 1ST  LAku  t l t>  FUkHATJ 


LOLGMNb  1-10  P AGE  kt  C GIVEN  E ENG  I k Ut  PLUItlNG  PAPtk  t Uk  1 PLUl» 

l.t.Y  ALUNL.  Ikt  HULL  ut  T kt  PAPtk, 

11-20  PAPtk*  s & 1 V t N * 1 0 I k Ut  PLUII1NU  PAPtk  t Uk  1 PLul, 
l.t.Y  ACkUSS  Ikt  kt.LL  it  Iht  PAPtk. 
r'l-iO  PSlAkl  s MINIMUH  value  Ut  THt  M-CLltikDlNA  (E  t uk  PluIUNG 
Ikt  k - A * I S . 

Jl-uu  lSIakT  s minimum  VALut  ut  Int  X-LUUHGlNAIt  tuk  PluIIIM. 
THt  i- A*Ib. 

«1-Su  Gt L Pli  s SCALE  PAkAMETE-kY  *kE  P E 1 I'»Ck  uf  PlLI  LEnGIm 

tGGAlS  GtLPU  INCHES  t.E  Ikt  VAHlAbLtS  k , / . 
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NtlLtP  PIStAPLH  INC 


! 


V tl»,  CM  If 


PAbt  MJ 


S 


• PAbAN  n A T A I NPU  1 LAPDS  (ONLY  Jf  iNPAbN  = 1) 

iHt  fULLUPlNb  bPUUP  Uf  CAHUS  PPOVlUtS  »N  UP1ION  f UP  lNPUlliNb  1 Ht  PAbAN 
tout  UUIPUI  DATA  InTu  POISSY  VIA  CAPuS.  Tbt  PAPTiCULAP  PtANlNb  Ut  THt 
UlfflPtNI  VAPlAhLlS  UlSLPlbtO  btLU"  IS  blVtN  IN  IHt  PAbAN  LUUt  USfcKS 
MANUAL.  ANU  SMUULU  ht  CONSUL  I tl1  PUP  iHllP  M f A n I N u # A SPItf  UtSCKiPUON 
* I LI  HI  b I V t N htkt  fUP  ClJNVItNLt, 


* PAbAN  HIAOIM.  CAPU 


f UPMA  T c 1 H A At  ) 


COLUMNS  1 - 7 1-  MtAU  = 7P  ALPMA/NUMtPlC  CHaPaCUPS  Uf  lIILt  IMUPMaTJOn 


• PAbAN  CUNIPUL  CAPU 


fUPMATlPlS.riO.S.Pflu.p.MO.S.SIS) 


CUL  UMNS 


l-S 

Nf  ILL 

S 

b-  1 U 

Nl 

2 

1 l-2u 

1 I Mt 

2 

21-30 

ALUNf 

2 

il-UU 

UIU 

2 

P 1 -SO 

SIPtCt 

2 

SI-SS 

K SUPf 

2 

Sb*bU 

N ft 

2 

b I "bS 

Apt 

2 

bb-70 

I M A X 1 

2 

7I-7S 

JMAX  1 

2 

f ILt  NUMtltP  Ut  UA  TA 
T I Ml  SItP  NUMBEP  Ut  (;  A 1 A 
ElASPStU  I I M t IN  IPAJttlUPY 
A L T I T UUt  IN  IPAJtdUPY 
NOSt  UP  vtLULIU 

CUPPtNT  S'AbNATIUN  P 0 1 N I PtCLLSIUN 
NUMBfcP  Ut  U b T SlUt  SUPtALt  PulNlS 
TUlAL  NUMhtP  UF  CUOKOlNAtt  API;  ItMP  PulNlS 
NUMhtP  Uf  bACKSlUfc  SUPtACt  PUINIS 
MAXIMUM  axial  CUUPUINATt  SPACiNb  INU1LY 
maximum  PAUIAL  COOPPlNAlt  SPACiNb  iNUicr 


• PAbAN  CtIUPt)  I NA  T t SPACiNb  UATA 


PUpma  USf  1 m . 1 1 J 


I Ht  PAUIAL  ANU  Axial  bASK  CUOPUlNAYt  SPACiNb  UaIa  A«t  G I v t N . S 1U  A CAPU, 
MlP  IPSPACt  l J ) , JS 1 , JMAX 1 ) ANU  l 7SPACI  ( 1 I , 1 = 1 , IMAX I ) ACLUPUINb  IU  Ibt 
Abuvt  POPMAI  UN  TWO  StPtPAIt  StPItb  Uf  CakuS,  Any  CaPu  NU I tlLCtU  IS  LtfT 
OLANX  flip  tACM  Stwits.  f UP  tXAMPLt,  f UP  Tht  PSPACt(J)  StPItb. 


COLUMNS 

1-lb 

PSPALt l 1 ) 

s 1 SI 

PAUIAL 

CUUPO 1 N A 1 1 

SM  A L 1 Mi 

DA  1 A 

1 7-32 

PSP»Lt (PJ 

2 

PAUIAL 

CUUPU I NA  Tt 

S^AC 1 NO 

U A 1 A 

33-UH 

pP-bU 

PSPALt (3) 

= 3KU 
2 ANU 

PAUI AL 

SU  UN  , 

LUUPU INA I L 

• • • • • 

SHAG INb 

DA  1 A 
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'tlLt*  KPSEAM.M  |NL..  "I,  VlC*.  lALlt- 


PAGAN  Pa  I a INPul  C»K|)3  CCu<!lMJtO) 


(UNL  1 If  1 NPAGN  S 1) 


* PAGAN  UU  t S 1 1-  fc  SGNPACE  Pull. I r.AKoS 


h.ki*«  I c 15. St  15,5) 


iMt  Uirsiu  SUHPACt  FUlMS  C LlflKO  I NA  1 1 S < TEMPtNAlUNE  A N L>  PKpSSUNE  OAlA  A«t 
blHNi  ONE  Tu  A CAKl)  f U»  ( ASUkP  ) PUJN1S  ACCUNOlNG  1U  THfc  AdUXE  PUKMAI, 


CliLuNNb  1-5  1 s SlmtaCt  NCUAC  POINT  NUMMEk 

O-P0  /*  s AXIAL  LUC'WG  I NA  I t OF  POINT 

PI-15  N*  = KALlAL  COUKUINaTE  Of  PUl  iv  T 

Jm-5U  1*  : I E P’Pfc  N A T UK  t up  PulNl 

51-05  PNS  = nuNkaL  PpESSUKE  ALIInG  A I POINT 

Oo-dU  PIS  z UM.tMUL  PNfcSSuWE  AC1IN G A T PlJlNl 


* pagan  I .y  rp  KJ  Ijk  ang  dACAOirt  SDWPALfc  POINT  LAHgS 


Pukka  1 ( J5.3P 15.51 


IMP.  I N 1 1 K 1 1 J W 4 0.1,  oACNSIDE  5IJK  F A C t POINTS  CUlIKO  I N A ! fc  S AND  TEmPEKATUKE  0 A 1 A 
Apt  GIVE*.  (JNt  Tu  A CANO  FUN  ( I = *SO KpPli,.,,P«  1 ACCUNOlNG  10  TmE  A0uvfc 
POKMAl,  KHtPp  THE  UUrtPtNlNl.  CONTINUES  PKUM  I M£  UUlSlllt  SuHPACE  PulNlS. 


columns  1-5  1 = iNltKluP  ON  b AC  *S  IDE  SuNPACt  Nul’AL  PoInI  NUMdEK 

o-PO  P*  r AXIAL  CUUHOINATE  OP  POINT 
pl-15  k-  s KAOlAL  COUNT)  Ina  ft  OF  PijInT 

Jo-50  Ik  = lEMPENA  T UP£  uF  POINT 


NO  I P . I«t  P A L.  A N NUUAl  PCJ 1 N 1 UAIA  IS  G l v t « dASIC*OLY  n Y 1 ME  LASI  I «U  SEKltS 
OF  CAKUS  In  UNt  LUNG  S1N1NG,  IN  1 Ht  UNUOK  UF  Cl)  Uulslut  SUNFALt  PulNlS, 
CP)  IMfPIUW  PulNlS  AND  (J)  KACNSIl/F  SUKPALE  Pi.iInTS,  ImE  PANIICULAK 
if  A P*  1 Nu  ANU  V 1 S IN  I MU  T IUN  AKt  FuCLY  fcXPLAlNtl)  IN  | ME  PAGAN  USt«S  MANuAI  . 


rEILER  RESEARlH  INC,,  Hi,  VILA,  CALIF 


page  no,  e 


XXXXXXXXXXXXXXXXXX 
POESSY  X COMMENT  PAGE  X 

xxxxxxxxxxxxxxxxxx 


RtlLER  RESEARCH  1 nC  , SULIL1TS  YOUR  COMMENIS  ABOUT  1 H 1 5 MANUAL  « l 1 M A VlER  lu 
IMPROVING  US  USEE  ULNLSS  IN  tATtR  EDITIONS, 


OU  YOU  E 1 NO  THIS  MANUAL  AOEUUATE  E UK  YUUEi  PURPOSE 


mha t improvement s du  you  recommend  to  better  serve  yuur  purpose 


NOTE  SPEC  I PIC  E PROPS  DISCOVERED  ( PLEASE  INCLOuE  PAUE  NUHbER  REFERENCE) 


GENERAL  COMMENTS 


FROM  NAME POSITIUN 

BUSINESS 

ADDRESS 


SEND  TO  rEILER  RESEARCH,  InC, 

SUITE  S2<4 

<<b72  BAYSHORE  FRONTAGE  RD. 
Ml,  VlER,  CALIFORNIA  R«uui 
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APPENDIX  A 


SAMPLE  PROBLEMS 

SAMPLE  PROBLEM  NO.  1 (PLUG  NOSE  TIP) 

The  first  sample  problem  choosen  to  illustrate  the  use  of  the 
POESSY  program  is  a typical  plug  nose  tip.  This  particular  plug  nose 
tip  was  analyzed  over  the  complete  trajectory  by  the  PAGAN  computer 
program  (Reference  4 ).  The  particular  point  choosen  to  illustrate 
the  use  of  the  POESSY  program  was  for  time  = 5.689  seconds,  altitude  = 
57,600  feet.  The  PAGAN  code  mesh  for  this  point  in  the  trajectory  is 
shown  in  Figure  6 . 

The  output  from  the  PAGAN  code  was  saved  on  a save  tape  and  thus 
was  used  for  input  to  POESSY,  i.e.,  TAPE  2 in  Figure  9 . Hence,  only 
three  cards  were  needed  for  input  to  POESSY.  A complete  description  of 
the  input  data  is  given  in  the  POESSY  User’s  Manual,  Section  VIII  of  this 
manual.  The  input  data  cards  for  this  plug  nose  tip  sample  problem  are 
shown  in  Figure  12,  where  (Card  1)  = Control  Card  and  (Cards  2 and  3)  = 
Plot  Scale  Control  Cards. 

One  will  notice  that  neither  DISIZ  or  DJSIZ  (Card  1,  Columns  46-65) 
were  specified  but  NELEM  (Card  1,  Columns  36-40)  was  specified  = 175. 
Hence,  the  sizing  of  the  basic  rectangular  element  was  performed  auto- 
matically. The  POESSY  program  printed  output  is  shown  in  Figure 
13.  Here  one  can  see  that  the  automatic  sizing  of  the  basic  rectangular 
element  resulted  in  DI  (=DR)  = 0.230  and  DJ  (=DZ)  = 0.230.  This  resulted 
in  producing  a final  structural  finite  element  mesh  having  223  nodal 
points,  192  elements,  17  outside  surface  pressure  boundary  conditions 
and  6 backside  reacting  surface  boundary  conditions.  The  resultant 
structural  finite  element  mesh  is  shown  plotted  (automatically  by  POESSY) 
in  Figure  14,  and  the  isotherm  contour  plot  in  Figure  15.  The  plot 
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scaling  information  for  these  plots  is  shown  in  the  printed  output  in 
Figure  13,  where  the  first  set  of  scaling  information  applies  to  the 
mesh  plot,  the  second  to  the  isotherm  contour  plot.  The  little  boxes 
shown  plotted  along  the  outside  boundary  of  the  mesh  in  Figure  14 
represent  the  intersections  of  the  PAGAN  mesh  with  this  boundary. 
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BAD  UUAO  REJECTIUN  PARAMETER  (THIFCT)  = -1S5.000 

dUUNUART  COnTUUR  TYPE  PARAMETER  UivdTYP)  a 1 

ORIGINAL  contour  PLOT  PARAMETER  (IORGCN)  - 1 

ANNOTaIION  parameter  (IANNQI)  : -E 

NUMBER  Of  CONTOURS  PARAMETER  (NCONT)  a 16 

MAXIMUM  PAPER  LEnGIH  (PAPERlJ  a 10,000 

MAXIMUM  PAPER  WIOIH  IPAPtRn)  A 11,000 

SCALING  PARAMETER  USCALt)  = 0 

ROTaIION  PARAMETER  (IRU1AT)  = 0 

PLOT  Margin  PARAME1ER  (IMA'rGn)  A 1 

minimum  COoRuINATE  value  parameter  (IMXMN)  = 0 

MINIMUM  R COORDINATE  VALUE  (RSTAR1)  = 0,000000 

MINIMUM  Z COORDINATE  VALUE  USTAR1)  a O.OuOOOO 

USER  SUPPLIED  SCALE  FACTOR  COELPO)  a 1.000000 


Figure  13.  POESSY  Output  for  Sample  Problem  No.  1 
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Figure  13.  Continued 
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Figure  13.  Continued. 
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Figure  13.  Continued 
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Figure  13.  Concluded. 
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SAMPLE  PROBLEM  MO.  2 (SHELL  NOSE  TIP) 

The  second  sample  problem  choosen  to  illustrate  the  use  of  the 
POESSY  program  is  a typical  shell  nose  tip.  Again  the  input  data  for 
different  points  along  the  trajectory  were  saved  on  a save  tape  and 
thus  were  used  for  input  to  POESSY,  i.e.,  TAPE  2.  Again  only  three 
cards  were  needed  for  input  to  POESSY. 

The  input  data  cards  for  this  shell  nose  tip  sample  problem  are 
shown  in  Figure  12,  where  again  (Card  1)  = Control  Card  and  (Cards 
2 and  3)  = Plot  Scale  Control  Cards.  One  will  notice  that  both  DISIZ 
and  DJSIZ  (Card  1,  Columns  46-65)  were  specified  as  DISIZ  = 0.150  and 
DJSIZ  = 0.275.  The  POESSY  program  printed  output  is  shown  in  Figure 
16.  Here  we  can  see  how  the  specified  basic  rectangular  element  size 
resulted  in  producing  a final  structural  finite  element  mesh  having  127 
nodal  points,  104  elements,  27  outside  surface  pressure  boundary  conditions 
and  4 backside  reacting  surface  boundary  conditions,  where  the  reacting 
pressure  acts  on  the  butt  end  of  the  conical  portion  of  the  shell  since 
the  parameter  ISHPLG  = 2 (cf.  Figure  12b,  Card  1,  Columns  6-10). 

The  resultant  plot  of  the  finite  element  mesh  is  shown  in  Figure  17 
and  the  isotherm  contour  plot  in  Figure  18.  Again,  the  plot  scaling 
information  for  these  two  plots  is  shown  in  the  printed  output,  Figure 
16. 


To  illustrate  how  the  same  shape  would  look  with  a more  dense  array 
of  elements,  this  same  sample  shell  nose  tip  problem  was  resolved,  using 
DISIZ  = 0.100  and  DJSIZ  = 0.124.  This  resulted  in  372  nodal  points, 

327  elements  and  59  pressure  boundary  conditions.  Only  the  plot  of  the 
resultant  finite  element  mesh  and  isotherm  contour  plot  are  shown  here 
in  Figures  19  and  20,  respectively. 
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CASE  NUMbE*  INUCASt)  a Ut 

S"ECL/PLUU  PARAMETER  USmPLu)  a i 

HECUESTEO  !!"£  STEP  NUMBER  (N5TEP)  a 126 
PRINT  QUT  PARAMETER  (JkkJI)  a 1 

PUNCm  CARD  P AR ame  T £ R UPunCm)  = 0 

PLLit  PARAMETER  ( I P UU  l ) = 3 

S A V L RESULTS  PAR  ARE  1 E R (IUSAVt)  A 0 

TOTAL  ELEMENTS  In  "ESN  Paha*eTeR  (nELEm)  - ISO 
r.AlA  INPOl  HARAKETtR  (lNPAGN)  A 0 

I-DIRECUUN  SPACING  PA"A"t!t<<  (ulSU)  = ,1500000 

J*D1RECTIUN  SPACING  Parameter  tUJSIZ)  - .2750000 

OAG  I.UAU  "EJEcUCN  parameter  (7RJFCT)  a *1  35,000 

BOUNDARY  CliNTuuR  7>Pt  PARAMETER  UMMYP)  A l 

UWIGINAC  CUnIUuR  Pool  PAHA"tTE«  UCRGCn)  I 1 

ANNOTATION  PARAMETER  (lANNCI)  a *2 

NUMbER  OP  CONTOURS  PARAMETER  (NCUM  1 a 16 

MAXIMUM  PAPtR  length  (PAPt"L)  a IN.OOO 

MAXIMUM  PAPER  "IGTh  IPaPEP")  s 11,000 

scaling  PaPANETEP  ( 1 SC  ALE ) * 0 

ROTATION  PARap.ETER  (IRCTaI)  a 0 

Plot  margin  parameter  UmarGn)  * 1 

MINIMUM  COORDINATE  VALUE  PARAMETER  (JMXMI,)  A 0 

MINIMUM  R CUURQ1NATE  value  (RSTaRT)  a 0,000000 

MINIMUM  z CUOROlNATE  value  ( ZST ART ) A 0,000000 

USER  SUPRLIEU  scale  FACTOR  (DLCPUl  a l.OoOOOo 


Figure  16. 
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Figure  16.  Continued 
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Figure  16.  Continued. 
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the  follu*.ing  is  the  nodal  point  information  fur  the  generated  mesh 
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Figure  16.  Continued. 
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Figure  16.  Continued. 
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Figure  16.  Continued. 
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Coarse  Mesh  Isotherm  Plot  of  Sample  Problem  No.  2 
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Figure  20. 


Refined  Mesh  Isotherm  Plot  of  Sample  Problem  No.  2 


APPENDIX  B 


PLOTTING  SOFTWARE 


There  are  two  types  of  plotting  software  referenced  in  POESSY, 

(1)  Subroutines  supplied  with  the  programs  pertaining  only  to  plotting 
functions  and  (2)  calls  to  what  may  be  called  a "standard  plotting 
software  package."  At  present  in  the  U.S.A.  there  are  no  established 
standards  for  plotting  software  as  there  are  for  language  compilers 
such  as  FORTRAN  IV,  COBAL,  etc.  Therefore,  in  the  author's  opinion, 
most  computer  installations  have  adopted  plotting  software  similar  to 
or  that  furnished  by  the  California  Computer  Products  Software 
(Reference  5)  commonly  referred  to  as  CALCOMP  software.  Calls  to  this 
software  is  what  is  incorporated  in  POESSY.  The  nature  and  type  of 
routines  contained  in  the  CALCOMP  software  will  be  briefly  explained 
here.  For  further  detailed  knowledge  of  these  routines,  the  reader 
is  referred  to  Reference  5. 

The  CALCOMP  Basic  Software  which  comprises  six  subroutines  are 
briefly  described  as: 

PLOT  converts  all  pen  movement  specification  from  inches 

to  actual  plotter  commands,  and  outputs  these  to  the 
attached  device. 

SYMBOL  draws  any  sequence  of  alphameric  characters,  as  well 

as  special  point-plotting  symbols.  It  calls  only  PLOT. 

NUMBER  draws  the  decimal  equivalent  of  an  internal  floating- 

point number.  It  calls  only  SYMBOL,  which  then  calls 
PLOT. 

SCALE  examines  a data  array  to  determine  an  optimum  starting 

value  and  a scalling  factor  for  use  by  AXIS  and  LINE  in 
converting  data  units  to  plotter  page  dimensions. 

It  is  the  only  subroutine  that  does  not  call  any  other, 
because  it  does  no  actual  plotting. 
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*1 


AXIS 


draws  an  axis  line  with  the  appropriate  scale  annotation 
and  title.  It  calls  SYMBOL  and  NUMBER  as  well  as  PLOT. 


L . 


LINE  plots  a series  of  scalled  data  points  defined  by  two 

arrays  (X  and  Y) , connecting  the  points  with  straight 
lines  if  desired.  It  may  call  SYMBOL  as  well  as  PLOT. 

and  the  actual  FORTRAN  calls  to  these  routines  is  given  by 
CALL  PLOT (XPAGE, YPAGE, I PEN) 

CALL  PLOTS (BUFFER, NBUF, ITAPE)  ***  ENTRY  POINT  IN  PLOT  ( ) 

CALL  FACTOR (FACT)  ***  ENTRY  POINT  IN  PLOT  ( ) 

CALL  WHERE (XP AGE, YPAGE, FACT)  ***  ENTRY  POINT  IN  PLOT  ( ) 

CALL  SYMBOL  (XPAGE, YPAGE, HEIGHT, IBCD .ANGLE ,NCHAR) 

CALL  SYMBOL (XPAGE , YPAGE , HEIGHT , I SYMNO , ANGLE , - ICODE ) 

CALL  NUMBER  (XPAGE, YPAGE, HEIGHT , FPNO .ANGLE, NDEC ) 

CALL  SCALE  (ARRAY,  AXLENG.NPTS, INC) 

CALL  AXIS  (XPAGE, YPAGE, I BCD, NCHAR , AXLENG , ANGLE, FI RSTV.DELTAV) 

CALL  LINE  (XARRAY , YARRAY ,NPTS , INC , LINTYP , I SYMNO) 

Again,  if  the  user's  plotting  software  differs  from  that  illustrated 
here,  he  is  referred  to  Reference  5 to  fully  understand  the  meaning  of 
all  of  the  arguments  in  the  calls  to  the  CALCOMP  software.  The  user  can 
then  make  the  appropriate  changes  in  POESSY  so  that  it  will  be  compatible 
with  his  plotting  software. 

The  subroutines  supplied  with  POESSY  that  pertain  strictly  to 
operations  associated  with  plotting  are  described  by 

PLTST  initializes  and  terminates  respectively  the  procedure 

PLTFN  by  calling  appropriate  CALCOMP  routines. 


**  * 
*★* 
*** 
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DASHLN 


TITPLT 


NPLABL 
ELLABL 
(NPLOC ) 


FRAMIT 


INNUM 

FPNUM 

(KPLACE) 

(STOREA) 


plots  a series  of  scalled  data  points  defined  by  two 
arrays  (X  and  Y) , connecting  the  points  with  a continuous 
dashed  line  of  length  specified  by  the  user  (DASH) . 

It  may  call  SYMBOL  as  well  as  PLOT. 

plots  a series  of  characters  (in  A4  format)  comprising 
a title,  given  the  coordinates  of  the  start  of  the 
title  and  appropriate  orientation  and  spacing  parameters. 
It  calls  only  SYMBOL  which  then  calls  PLOT. 

annotates  a boundary  contour  line  defined  by  a series 
of  contiguous  nodal  points  (NPLABL)  or  a series  of 
element  sides  (ELLABL)  with  the  appropriate  nodal 
point  or  element  number,  placed  next  to  the  point 
or  line  they  are  describing.  They  call  auxiliarly 
routines  NPLOC  and  INNUM  as  well  as  SYMBOL,  WHERE 
(PLOT  entry)  and  PLOT. 

determines  all  of  the  required  sizing  and  scalling 
data  necessary  to  fit  a plot  onto  a specified  paper 
size  (length  and  width) , given  the  parameters  controlling 
the  plot,  such  as,  title  information,  margin  specifi- 
cations, rotation  and  scaling  specifications  and 
maximum  and  minimum  coordinate  values  of  the  points 
defining  the  plot.  It  calls  an  auxiliary  routine 
PSCALE,  which  scales  both  axis  of  the  plot  producing 
the  axes  lengths,  starting  values  and  a scale  factor 
for  use  by  AXIS. 

changes  an  integer  or  floating-point  number  from  internal 
computer  based  format  into  alpameric  format  in  the  forms 
Im  (INNUM)  and  Fm.n  (FPNUM)  following  standard  FORTRAN 
output  conventions.  They  call  auxiliary  routines  KPLACE 
which  locates  the  word  and  position  in  the  resulting 
alphameric  character  string  and  STOREA  which  stores  a 
character  in  a computer  word. 


The  type  of  annotation  (numbering)  performed  by  the  subroutines 
NPLABL  and  ELLABL  is  shown  graphically  in  Figure  21.  One  will  notice 
that  the  element  numbering  is  placed  at  the  mid-point  along  the 
element's  side.  Also,  the  nodal  point  numbering  is  orientated  by 
the  average  of  the  outward  normals  of  the  two  connecting  sides. 

The  dashed  lines  in  Figure  21  are  only  included  to  illustrate  the 
spacing  and  are  not  drawn  by  the  subroutines  NPLABL  and  ELLABL. 

The  rest  of  the  subroutines  are  self-explanatory,  with  sufficient 
comment  cards  included  in  the  coding  to  explain  the  meaning  of  all 
operations. 
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If  a user  does  not  have  access  to  the  Basic  CALCOMP  Software, 
he  may  (a)  easily  change  all  of  the  calls  to  this  software  to  appropriate 
calls  to  his  own  plotting  software  which  performs  the  same  functions  or 
(b)  simulate  the  CALCOMP  routines  with  dummy  subroutines  which  then  call 
his  appropriate  software.  The  author  apologizes  for  the  fact  that  no 
standards  exist  for  plotting  software,  however  he  feels  that  the  CALCOMP 
Basic  Software  comes  the  closest  to  all  those  being  used  at  the  present 


time. 


Boundary  Contour  Line 


Character  Height  = 0.085  in. 
DOUT  Spacing  = 0.210  in. 


Figure  21.  NPLABL  and  ELLABL  Annotation 
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